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ABSTRACT 

We shall present results of a recent collabora- 
tive effort between the authors attempting to im- 
plement the numerical flux scheme, AUSM + and 
its new developments, into a widely used NASA 
code, OVERFLOW. This paper is intended to 
give a thorough and systematic documentation 
about the solutions of default test cases using the 
AUSM+scheme. Hence we will address various 
aspects of numerical solutions, such as accuracy, 
convergence rate, and effects of turbulence mod- 
els, over a variety of geometries, speed regimes. 

We will briefly describe the numerical schemes 
employed in the calculations, including the ca- 
pability of solving for low-speed flows and mul- 
tiphase flows by employing the concept of nu- 
merical speed of sound. As a bonus, this low 
Mach number formulations also enhances con- 
vergence to steady solutions for flows even at 
transonic speed. Calculations for complex 3D 
turbulent flows were performed with several tur- 
bulence models and the results display excellent 
agreements with measured data. 

1. INTRODUCTION 

A collaborative attempt has been evolved to im- 
plement the numerical flux scheme, AUSM + and 
its new developments [1-5], into a NASA CFD 
code, OVERFLOW [6]. The OVERFLOW code 
currently is a production code that has been 
widely used in the US industry, government labs, 
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and universities. It is written to use the over- 
set, structured (Chimera) grid systems and has 
evolved through several significant changes, due 
to contributions from various researchers over 
more than two decades (see [6] for further details 
and references therein). The Chimera grid ap- 
proach was developed for multibodies and com- 
plex geometry flow simulations. Major applica- 
tions have included the Space Shuttle Launch 
Vehicle and subsonic transport high-lift config- 
urations. The Pulliam-Chaussee diagonalized 
implicit algorithm [16] has become the stan- 
dard solver for the code. Several numerical flux 
schemes are commonly used, including central 
differencing and Roe upwind differencing. 

The AUSM family schemes, even though rela- 
tively young, also has gone through several mod- 
ifications. The schemes have been made into 
some major production codes as well (see, for 
example, [7]-[9]). 

In addition, this collaboration is intended to give 
a thorough and systematic documentation about 
the solutions of default test cases obtained by us- 
ing the newly implemented scheme. Hence, it is 
considered to be a long-term project, especially 
the validation processes. As a result, this paper 
contains only the first installment of results from 
this continuing effort. The numerical schemes 
will be briefly described, together with the con- 
cept of numerical speed of sound [5], which is 
especially useful for low speed flows and multi- 
phase flows [11,17]. To confirm the effectiveness 
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of the AUSM + scheme for practical engineering 
problems, calculations for a wide variety of flows 
were made. The results display excellent agree- 
ments with measured data. 

2. NUMERICAL METHODS 


Implicit Operator 

The nonlinear equations of (4) are linearized, 
resulting in an implicit equation to solve for 

AQ = Q n+1 - Q n . 


Governing Equations 

We solve the Reynolds-averaged Navier-Stokes 
equations written in curvilinear coordinates 
(£, and conservation form: 

Qt + Fu + F^ = 0, 1= 1,2,3 (£, 17 , 0 , (1) 
where 

Q = (p,pu,pv,pw,pe t ) T /J. 

The inviscid and viscous fluxes are included in 
F/ t / and f{^, where the summation and differ- 
entiation notation, F ij = Ya=\ dFi/dxi is as- 
sumed. 
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The flow variables in Q and the geometric Ja- 
cobian ./ are standard. The turbulence terms 
are included in F^ and expressed in the form 
of eddy- viscosity described by algebraic or dif- 
ferential equations. 

To enhance the convergence rate when solving 
flows at low speeds, the time-derivative term is 
premultiplied by a conditioning matrix T. The 
equations now can be cast as, 

TQ ( + F /t , + FjV = 0, Q = (p, «, v, w,T) t /J. 

(4) 

Several forms for the local preconditioning ma- 
trix r have been proposed, notably by Choi and 
Merkle [12], Turkel [13], and Van Leer, Lee, and 
Roe [14]. While these matrices are different, ow- 
ing to different design criteria, it is not essential 
for now to explicitly specify the form of T. 
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Note that the LHS implicit operator consists of 
only the inviscid part of the whole flux. Several 
methods have been suggested for approximating 
the implicit operator, in order to gain conver- 
gence rate, reduce computation, and increase 
stability. Several options are available in the 
OVERFLOW code. At the initial stage of this 
study, our focus is mostly on the accuracy and 
convergence rate of steady solutions, which are 
controlled only by the numerical methods repre- 
senting the RHS and should be independent of 
the LHS operator because it vanishes as steady 
solutions are reached. In this study, we fixed 
the implicit operator to be the Beam Warming 
ADI scheme [15], in which the LHS was approx- 
imately factored according to the computation 
coordinates in the following form: 


[I + Affix NAj] [I + Atf 2 NA 2 ] [I + At.fi 3 NA 3 ] AQ 
= -AtNS, (f, + f{ v) ) N = Mf'. 

( 7 ) 

The LHS operator was then discretized using 
the central differencing method, with appro- 
priate second-order and fourth-order smooth- 
ing terms added for stability. Furthermore, the 
diagonalization procedure of Pulliam-Chaussee 
[16] was applied. Details can be found in [6]. 
Since the RHS operator was represented using 
the AUSM + scheme - an upwind scheme, there 
should be a concern about the inconsistency in 
discretization of the LHS and RHS terms. And 
the stability and convergence rate are likely not 
at optimal situations. This, however, is a very 
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difficult issue by itself and will not be addressed 
in the present paper. 

Explicit Operator 


In this section, we will give a brief description 
of the algorithm involved in the AUSM + [2] and 
the new development for the low-speed flows 
[5,11,17]. Concerning detailed numerical prop- 
erties and detailed analysis, the reader should 
consult the cited references. 

For illustration purposes, we shall begin by con- 
sidering the ID flux for ideal gas. The inviscid 
flux F is writ ten as a sum of convective and pres- 
sure fluxes: 


= pu 



: rn 4>+p. 


(8) 

The vectors <3> and p are introduced respectively 
as shown. It is noted that a common mass flux 
m appears in all equations. This is also true for 
multidimensions. Since the mass flux is common 
for all equations, its effects will thus perpetuate 
in all variables. Hence, it is desirable observing 
this fact at the discrete level when devising a 
new scheme. However, this fact is not entirely 
enforced in several modern numerical schemes. 
If this same factor is kept at the discrete level 
for all equations, it becomes easy to add other 
conservation laws. 


As observed in (8), all one needs to do is to define 
the quantities, ra, <I> and p. Concerning the de- 
tails of developments of the AUSM-family fluxes, 
the interested readers are referred to the refer- 
ences [1-5,11]. In what follows we will summarize 
only the latest formulas implemented in the code 
[5]. Firstly, a simple upwinding is applied to the 
convected quantity to yield the numerical flux 
at the cell interface (denoted by subscript 1/2) 
straddling cells j and j + 1. 


fl/2 
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The pressure flux is nothing but 


PI/2 = + V~ a) (M ]+ x ) Pj+l , ( 10) 
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The only task left now is to find the mass flux 
m, which can be easily done by following these 
steps. 


1. A concept of numerical speed of sound is 
introduced in [5] to make the flux formu- 
lation valid for the entire speed regime. 
Here, we let the speed of sound be scaled 
by a scaling factor f a (M ; A/*), 


j2 — fa ( A/ , A/* ) A 1 /2 1 

(12) 

where 
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A/ = (Mj + Mj+\ )/2, 

(14) 

and 


A/ 2 = min(l,max(A/ 2 , A/ 2 0 )). 

(15) 

The unsealed speed of sound ct \/2 can 

take 


several forms, as given in [2,4], one of 
which allows exact capturing of stationary 
shock. For simplicity, a simple average suf- 
fices, 

a | /2 = [dj + a J+x )/2. (16) 

The cutoff parameter M co << 1 is intro- 
duced to prevent a singularity at stagna- 
tion point. It is a user-specified parame- 
ter and we choose A/ 2 0 = A/^/4. The use 
of numerical speed of sound a results in a 
reduction in numerical dissipation at low 
Mach numbers. 


2. Use a to define 
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This scaled Mach number will revert to 
the local physical Mach number at super- 
sonic speeds. In fact, the scaling factor be- 
comes essential only in the low Mach range 
since it rises rather quickly towards unity 
as M increases. Hence, what the numerical 
speed of sound does is to make the Mach 
number appear larger in determining the 
numerical fluxes in the low speed regime. 

3. Let 

Mj = h(l + Ml)Mj + (1 - M, 2 )M j+1 ], 

( 18 ) 

M J+ i = h(l + M*)Mj+i + (1 - Ml)M 3 ). 

(19) 

4. Set the interface Mach number, 

M 1/2 = M+ 3) (Mj) + 1 )- 

(20) 

The split Mach number functions, 3 y 
are defined as follows: 

' i(M±|M|), if |A/| > 1, 
Mf 4j)) (M) = | ±\(M ±l) 2 ± §(M 2 - l) 2 ; 

otherwise. 

( 21 ) 

The numeral in the subscript of ^ t 
indicates the degree of polynomials used. 

5. The mass flux can now be obtained by a 
simple upwind formula, 

"*1/2= + P}+ 1 ) 

+ |^l/2l (Pj “ ^ + 1)]’ (22) 


at once, the moving shock around a 90-degree 
corner is a suitable problem. Several modern 
upwind methods have been known to encounter 
difficulties [4]. Figure 1 displays the density con- 
tours, showing crisp resolution of the original 
and several internal shocks, as well as the shear 
layers resulting from the intersection of shocks. 

itvd= 2 muscl: primitive variables 



Fig. 1 Supersonic corner problem. Calculation 
was made with linear extrapolation of primitive 
variables using the Van Albada limiter 

To see the effect of this numerical speed of sound, 
we evaluate the condition number at low speed, 

b= M+“_0(1), ,23) 

\u\ 


For identification purpose, we now call this ex- 
tended method, AUSM + -a, to highlight the role 
of the numerical speed of sound a. 

Remark : Test cases have shown that the defi- 
nitions in Eqs. ( 18 )-( 19) are not strictly nec- 
essary, the simpler alternatives Mj = Mj and 
A/j+i = Mj+\ can be used as well. 

To demonstrate the robustness, shock-capturing 
capability, and positivity-preserving property all 

tThe coefficient in (19b) of 7C/T29, 364-382 (1996), s] 


because 

a -+ 0, as 1 1 / 1 — > 0. (24) 

Hence, the condition number remains of order 
unity for all speed ranges. Also, the numerical 
dissipation based on this new speed of sound now 
scales with local speed |u|, instead of local speed 
of sound a as |u| — 0. 

ould be 1/4, not 1/2. 
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To further improve the residual convergence in 
the low speed range, it is found beneficial in Ed- 
wards and Liou [11] if a pressure diffusion term 
m p is included in the mass flux. 

m 1 / 2 = Eq. (22)+ m p . (25) 


Here rh p is expressed in a general form, 

'»p= 1 - M* )&M (Pj - Pj+i ), (26) 

where 


A M = 






(Mj) - MfoiMj) 

) + '+,)(Mj+l )] ’ 


(27) 


and the function V in the denominator can take 
several forms. Based on the mass flux of the 
AUSMDV scheme, Edwards and Liou [11] de- 
rived the pressure diffusion term for low Mach 
nu mbers, 


V = (28) 

l Pj Pj+ i J 

Another formula, has been proposed [5], 

T> = — \M*(pj + p J+ 1 )(pj + pj . h ) 

PjPj+i L 

-(pj ~ Pj+\ )(Pj ~ Pj+i )] • (29) 


Furthermore, the last term can be omitted to 
guarantee that V be positive, 


V 


= — - — \ M i(vj + pj + 1 )(pj + Pj+ 1) 
pjpj+i 1 


(30) 


which looks similar to the first expression. The 
scheme with the inclusion of the pressure diffu- 
sion term (Eq. (22)) is now denoted as AUSM + - 
ap. 

We now r make remarks on the preconditioning 
matrix T. We have used the Weiss-Smith T [18] 
to arrive at the scaling function f a (M: A/*) in 
Eq. (12). Other preconditioners [12-14] can be 
used as well. The procedure for extension will 
be precisely the same since all one needs is the 
eigenvalues of the preconditioned hyperbolic sys- 
tem. Thus, a new numerical speed of sound a 


can be expressed in terms of the scaling function 
f a ( M ; A/* ). However, no significant effect on the 
solution is anticipated because all these precon- 
ditioners yield more or less the same behavior in 
the limits of M — ► 0 and 1. Unless at low speed 
(say Moo < 0.3), it was found in our calculations 
not necessary to include the preconditioning ma- 
trix in solving the governing equations. In other 
words, the scaling function can be incorporated 
alone, as in Eqs. ( 12 )-( 22 ), in the numerical flux 
and improvements in accuracy and convergence 
can be realized. 

3. RESULTS AND DISCUSSION 

In this section, we will present 2D and 3D 
Navier-Stokes solutions for turbulent flows over 
various geometries. The scheme proposed in 
this paper was implemented in the OVERFLOW 
code that has been developed by Bulling et al. 
[6]. The LHS operator was approximated with 
the standard central difference scheme plus ap- 
propriate artificial damping terms, (even though 
the RHS residual operator was represented with 
an upwind scheme !), it was then further fac- 
tored and diagonalized in each space dimension. 
To accelerate convergence, a full multigrid strat- 
egy using two levels of coarser grids is applied 
initially. 

The flux in the RHS operator was constructed 
with a third-order accurate interpolation for the 
primitive variables, together with limiter used 
by Koren [19]. The cutoff Mach number in Eq. 
(15) is given bv M^ 0 = A/^/4. Also, the OVER- 
FLOW code has a parameter controlling the use 
of the preconditioning matrix T. We kept the 
default value to be 3A/^ < 1 under which con- 
dition T was activated. 

For all calculations presented in this paper, the 
flow is assumed fully turbulent — no laminar or 
transition regions are consider. Unless stated 
otherwise, the baseline model is the Spalart- 
Allmaras one-equation model. But the effects 
of various turbulence models on the solutions 
in comparison with measured data will also be 
studied. 

In this paper we will demonstrate the effective- 
ness of using the numerical speed of sound in 



calculating flows at all speeds, specifically focus- 
ing on two issues: (1) convergence rate and (2) 
accuracy. 

NACA Airfoil 

A C-grid of 249 X 56 points was used for the 
NACA 0012 airfoil. The purpose of calculating 
this flow was to determine the performance of 
five popular turbulence models against the ex- 
perimental data, see [20], so that it may help 
down-select fewer models for additional vali- 
dation tests. These five models include the 
Baldwin- Lomax (denoted as BL) [21], Baldwin- 
Barth (BB) [22], Spalart-AUmaras (SA) [23], 
Mentor k-c (SST) [24], and k-u? [25] models. Fig- 
ure 2 displays the surface pressure coefficient on 
the airfoil, indicating a relatively strong shock 
situating in the middle of the suction surface of 
the airfoil. In fact the boundary layer is sep- 
arated because of the adverse pressure set up 
by the shock. The Baldwin- Lomax model un- 
derestimates the effects of shock/boundary layer 
interaction, resulting in an incorrect shock loca- 
tion; but it gives the best result on the pressure 
side. The Baldwin- Barth one-equation model 
gives a very good prediction of the shock loca- 
tion, but with some disagreement with the data 
on the pressure side. The Spalart-Allmaras one- 
equation model overall gives a good agreement 
with the data. The suction-side pressure is best 
predicted by the SST model, but the pressure 
side is accurately predicted by the BL and k- 
u) models. Figure 3 displays the pressure con- 
tours resulting from all five turbulence models. 
Noticing that the models that predict the shock 
location correctly appear to have a shock pro- 
file curving forward more pronouncedly at the 
foot. This behavior also corresponds to the fact 
of having a larger separated bubble. 

The calculations were carried out for all models 
using the same set of parameters associated with 
time marching procedures and the implicit oper- 
ator, which are by no means the optimal set for 
any models. The convergence histories of these 
five models are shown in Fig. 4, where the SA 
model reaches the steady value of Cl in about 
2,000 time steps, while the BL model takes about 


3,700 steps. It is interesting that the lift coeffi- 
cients from the BL and BB models starts out at a 
low value and increases to an asymptote as iter- 
ation proceeds and a completely opposite trend 
is observed for the other models. 

Shuttle External Tank 

This is an axisymmetric Shuttle external tank 
geometry with a sharp nose and blunt base, 
downstream of which there is a significant sep- 
aration zone. One of the grid lines conforms to 
the body and grows outward and a plane con- 
sists of 88 x 60 grid points. Shown also are the 
meshes clustered around some key regions, one of 
which is in the middle to resolve a tiny notch (not 
visible to the scale). The free stream Reynolds 
number was fixed at 10,000. We have tested con- 
ditions from low Mach, transonic, to supersonic 
flows. Several schemes were considered, consist- 
ing of the standard AUSM + , AUSM + -a, and 
AUSM + -ap, with and without the Weiss-Smith 
preconditioner. In all calculations for this prob- 
lem, we made 200 steps for each of two coarser 
grids prior to the finest grid, on which 3000 more 
steps were continued unless noted otherwise. 

Table I summarizes the convergence behavior of 
the above combinations. We observe the fol- 
lowing: (1) For low Mach numbers (approxi- 

mately A/co < 0.3), it was found necessary to use 
the time-derivative preconditioner T so that the 
numerical dissipations in both the implicit and 
explicit operators are compatibly scaled. Oth- 
erwise, the calculation either diverged or stag- 
nated. (2) For flows at transonic speeds or 
higher, the time-derivative preconditioner, as 
given in Eq. (4), serves no benefits whatso- 
ever, even though the fluid speed is low in the 
viscous and the base flow regions. This is un- 
derstandable because the preconditioner effects 
only the inviscid waves and the information in 
the viscous-dominated regions is only transmit- 
ted via diffusion processes which are ably han- 
dled by the implicit operator. 

In Fig. 5, we display the convergence history 
for various Mach numbers using AUSM + , but 
without the preconditioner T. The residuals for 
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the low Mach-number cases stall after a drop 
of four orders of magnitude. These drops in 
many calculations, although not especially ad- 
mirable, would have been acceptable. However, 
a close examination of the solution reveals that it 
is completely unacceptable, as shown in Fig. 6. 
It appears that there is a false boundary (exactly 
aligned with a grid line) at which information is 
unable to get passed. This phenomenon is quite 
typical in the low Mach-number calculations us- 
ing an unmodified compressible code, also seen 
in [5]. Hence, a measure of caution should be 
taken when reading the residual history for the 
low Mach-number solutions. 

On the other hand, the convergence histories 
with the use of the numerical speed of sound 
display improvement over those without it, as 
shown in Fig. 5. As noted earlier in Table I, it is 
necessary to invoke T for the Afoo — O.Oland 0.1 
cases. The convergence rates for these two cal- 
culations nearly coincide with each other, indi- 
cating Mach-number independence. 

In Fig. 7, we show the solution at N— 1000 steps 
at which the residual has been dropped to the 
level approximately equal to that of the baseline 
scheme AUSM+, at N=6400. However, the so- 
lution now is well behaved and is every bit as 
good as the final solution at N=3200 at which 
the residual has been further reduced by two or- 
ders. Also the blow-up view near the surface re- 
veals smooth profiles of pressure contours, unlike 
the standard AUSM + which has been known to 
yield unwanted pressure oscillations in viscous 
layers along the transverse grid lines when the 
mesh aspect ratio is large and flow is essentially 
parallel to a grid line.* 

Finally the effect of including the pressure diffu- 
sion term on the solution was investigated and 
the convergence histories are also included in 
Fig. 5. Again, the preconditioner T must be 
used for the low Mach-number cases and their 
convergence histories are essentially identical, 
becoming independent of Mach number as the 
Mach number lowers. The pressure contours are 
indistinguishable from those shown in Fig. 7 and 
are thus not included. 


Comparing results in Fig. 5, we see that the 
convergence rate is improved in the transonic 
ranges bv simply using the numerical speed of 
sound alone. For low Mach number cases, M ^ — 
0.01.0.1, another order of reduction can be ob- 
tained by including the pressure diffusion term. 
Also, the use of numerical speed of sound yields 
the convergence histories that are relatively in- 
sensitive to the flow speeds. 

We now summarize major findings from the 
study of this problem: (1) The numerical speed 
of sound concept is an effective means of ex- 
tending the AUSM-type discretization to solve 
low Mach number flows in an accurate and effi- 
cient manner. (2) Since the numerical speed of 
sound is reduced with the flow speed, the nu- 
merical dissipation changes accordingly, and a 
compatible implicit operator (one that includes 
the preconditioning matrix) must be used. (3) 
For moderate Mach numbers and beyond, it is 
not necessary to use T, (4) Incorporation of the 
numerical speed of sound, as described in steps 
(12)-(22)), helps remove pressure oscillations in 
the viscous layers. 

ONERA M6 Wing 

The next problem is the ONERA M6 wing with 
the free stream conditions M iyo = 0.84, and 
Re 0 o — 18.2 x 10 6 under various angles of attack. 
The computation domain consists of 269x35x67 
grid points. For this case, the precondition- 
ing matrix T was automatically turned off in 
the code since the controlling factor 3 A/jv ex- 
ceeds unity. However, the numerical speed of 
sound a was active with M co = A/^,/2. The 
pressure contours on the wing surfaces are dis- 
played in Fig. 8, showing the well-known A- 
shock pattern, where the double-shock region 
extends more than 80% of the wing span. The 
detailed comparison of surface pressure distribu- 
tions are shown for three spanwise sections, Figs. 
9, 10 and 11 respectively for 44%i, 65%;, and 90% , 
for three angles of attack. The computed results 
are in very good agreement with the data [26], 


* However, the pressure distribution along the wall is smooth. 



noticing the accurate capturing of the shock lo- 
cations. Recall that the SA model is assumed, 
unless noted otherwise. 

The convergence histories are presented in Fig. 
12 for three different angles of attack. They show 
a continuing decrease in 800 steps (600 fine-grid 
steps) by about four orders of magnitude from 
the largest values, at nearly the same rates even 
though the flow features have varied consider- 
ably. 

Figure 13 compares the effect of three turbulence 
models, SA, k — u;, and SST models, on the pres- 
sure prediction; the SA and SST models give 
nearly identical results, except minor differences 
behind the shock at rj = 0.95, and the k — 
model seems to give a slightly better agreement 
with the data for the shock location. The parti- 
cle traces of two angles of attack, a — 3.06 and 
5.06, at different wingspans are given in Fig. 14, 
showing little spanwise excursion, except at the 
wing tip for the a = 3.06° case, but appreciable 
three-dimensional effect for the other case. 

Wingbody 

Turbulent flows over a wingbody configuration 
were calculated. The geometry is shown in Fig. 
15, where the sting is included in the calculation. 
The computation domain was gridded using the 
chimera overset grid technique and entire grid 
composed of seven grids. The particle traces 
around the body and wing, shown in Fig. 16, 
indicate that the flow follows the configuration 
quite closely and is essentially two-dimensional 
over the wing. Figure 17 depicts the pressure 
coefficients at various spanwise locations. The 
computed results are in excellent agreement with 
the measured data [27]. Moreover, the pressure 
coefficients along the body, shown in Fig. 18, 
exhibits the similar level of excellent agreement 
with the data. Even in the wing root region 
where a sharp variation is encountered. 

Finally, Fig. 19 displays a well-behaved con- 
vergence history, reducing the residual error by 
more than five (5) orders of magnitude in 800 
steps. 

4. CONCLUDING REMARKS 


We have successfully implemented the new 
AUSM+flux scheme in the production code of 
NASA, OVERFLOW and presented in this pa- 
per the validation results of test cases. With the 
introduction of the numerical speed of sound in 
the previous AUSM+flux scheme, it now yields 
accurate results at low Mach number and effi- 
cient convergence at a rate (nearly) independent 
of Mach number and angle of attack. Inter- 
estingly, the convergence rate is even enhanced 
at transonic speeds with the numerical speed of 
sound, without invoking the preconditioning ma- 
trix. Effect of turbulence models on the results 
has also been investigated, the SA and SST mod- 
els give the best results in comparison with the 
measure data of NAGA 0012 airfoil and ONERA 
M6 wing. 
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Table 1: Summary of convergence behavior due to various schemes for the shuttle external tank, 
Re^ = 10000 . 


" Moo 

Scheme 

No Precond. 

Precond. 


AUSM+ 


Diverg. 

0.01 

AUSM+-a 

No converg. 

y/ 


AUSM + -a + m p 

No converg. 

x/ 


AUSM+ 

v/ 

V 

0.80 

AUSM+-a 

V 

y/ 


AUSM+-a + m p 

y/ 

V 


AUSM+ 

V 

V 

2.00 

AUSM+-a 

V 

V 


AUSM+a -(- nip 

V 

V 
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Fig. 5 Convergence history for the shuttle external tank problem: (Left) AUSM + , (Middle) ALSM + -a 
(Right)AUSM+-ap. 



Fig. 6 Pressure contours for the shuttle external tank problem obtained at N— 6400 time steps for — 0.01, 
using the standard AUSM + . The picture to the right shows a magnified view near the nose. 
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Fig. 7 Pressure contours for the shuttle external tank problem obtained at N— 1000 time steps for A/^c — 0.01, 
using AUSM + -a. The other picture shows a magnified view near the nose. 
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Fig. 8 Pressure contours on the pressure and suction surfaces of the ONERA M6 wing at a = 3.06°, A/oo = 
0.84, and Re c c = 18.2 x 10 b , showing the A-shock pattern near the wing tip. 



Fig. 9 Pressure distribution at the spanwise section, 44% of the ONERA M6 wing, for various angles of attack, 
A/co = 0.84, Reco = 18.2 x 10 6 . 



Fig. 10 Pressure distribution at the spanwise section, 65% of the ONERA M6 wing, for various angles of attack, 
A^c = 0.84,/Jcco = 18.2 x 10 6 . 



Fig. 1 1 Pressure distribution at the spanwise section, 90% of the ONER A M6 wing, for various angles of attack, 
Moo = 0.84, Reoo = 18.2 x 10 6 . 
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Fig. 12 Convergence history for the ONERA \16 wing problem at three angles of attack, M x = 0.84, Re^ = 
18.2 x 10 6 . 


a. 

U 



Fig. 13 Effect, of turbulence models on the prediction of pressure distributions for the ONERA \16 wing problem 
at a - 3.06°, M~o — 0.84, Re*, = 18.2 x 10 6 .; Solid line: SA model, dash line: model, and small dot: SST 

model. The SA and SST results are nearly identical, except some differences behind the shock at. ij = 0.95. 
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Fig. 18 Pressure distribution at angles along the body, A/oo — 0.8, a = 2 Re^ — 0.167 x 10 b . 





